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Abstract 

We  consider  the  analysis  of  linear  programming  relaxations  of  a  large  class  of  combinatorial 
problems  that  can  be  formulated  as  problems  of  covering  cuts,  including  the  Steiner  tree,  the 
traveling  salesman,  the  vehicle  routing,  the  matching,  the  T-join  and  the  survivable  network 
design  problem,  to  name  a  few.  We  prove  that  all  of  the  problems  in  the  class  satisfy  a  deep 
structural  property,  the  parsimonious  property,  generalizing  earlier  work  by  Goemans  and  Bert- 
simas  [3].  We  identify  two  set  of  conditions  for  the  parsimonious  property  to  hold  and  offer 
two  proof  techniques  based  on  combinatorial  and  algebraic  arguments.  We  examine  several  con- 
sequences of  the  parsimonious  property  in  proving  monotonicity  properties  of  LP  relaxations, 
giving  genuinely  simple  proofs  of  integraHty  of  polyhedra  in  this  class,  offering  a  unifying  under- 
standing of  results  in  disjoint  path  problems  and  in  the  approximabihty  of  problems  in  the  class. 
We  also  propose  a  new  proof  method  that  utilizes  the  parsimonious  property  for  establishing 
worst  case  bounds  between  the  gap  of  the  IP  and  LP  values.  Our  analysis  unifies  and  extends 
a  large  set  of  results  in  combinatorial  optimization. 
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1      Introduction 

We  consider  the  following  class  of  problems  defined  on  a  graph  G  =  {V,E)  and  described  by  the 
following  integer  programming  formulation 

IZf{D)  =    minimize  SeCE^ea^e 

subject  to    Ee€*(t)  ^e  =  f{i),  i  e  D  CV 

T,e€S(S)^e>f{S},SCV 

where  /  :  2^  -+  Z+  is  a  given  set  function,  6{S)  =  {e  =  {i,j)  e  E\  i  e  S,  j  €  V  \  S}.  By 
selecting  different  set  functions  f{S)  and  different  sets  D  we  can  model  a  large  class  of  combinatorial 
problems,  including  the  Steiner  tree,  the  traveling  salesman,  the  vehicle  routing,  the  matching.  T- 
join  and  survivable  network  design  problem,  to  name  a  few. 

Let  IPf{D)  be  the  underlying  feasible  space.  We  denote  the  LP  relaxation  as  Pf{D),  in  wliich 
we  replace  constraints  Xe  €  Z+  with  Xe  >  0.  We  denote  the  value  of  the  LP  relaxation  as  Zf{D). 
Goemans  and  Bertsimas  [3]  studied  the  survivable  network  design  problem,  in  which  the  objective 
is  to  design  a  network  at  minimal  cost  that  satisfies  connectivity  requirements  (for  each  pair  {i.j) 
of  nodes  in  V,  the  solution  should  contain  at  least  Tjj  edge  disjoint  paths)  and  considered  an  integer 
programming  formulation  of  the  type  IPj{D)  with  f{S)  =  max(,_j)gi(5)  rjj,  D  =  0.  They  showed 
the  following  property,  which  they  call  the  parsimonious  property,  which  in  our  notation  can  l)c 
stated: 

Theorem  1  [Goemans  and  Bertsimas  [3]]  //  the  costs  Ce  satisfy  the  triangle  inequality  {c,j  < 
Cifc  +  Ckj  for  all  i,j,  k  €  VJ,  then  for  the  survivable  network  design  problem  {f{S)  =  max^^^^^s)  ^e) 

Zj{D)  =  Zj{%). 

In  other  words,  the  degree  constraints  are  imnecessary  for  the  LP  relaxation  in  the  stirvivablc  net- 
work design  problem.  They  further  examine  several  sometimes  surprising  structural  and  algorithniic 
properties  of  the  LP  relaxation,  and  examine  the  worst  case  behavior  of  -^^4^  for  the  survivable 


network  design  problem.  Goemans  and  Williamson  [4],  Williamson  et.al.  [16]  and  Goemans  et.  al. 
[5]  show  interesting  worst  case  bounds  on  the  ratio  -^^my- 

Our  goal  in  this  paper  is  to  understand  the  class  of  problems  for  which  the  parsimonious  propertj' 
holds  and  examine  several  implications  of  the  parsimonious  property.  In  this  way  we  shed  new  light 
to  a  large  collection  of  results  in  discrete  optimization  and  graph  theory,  imderstand  their  common 
origin  and  generalize  them  in  interesting  ways.  In  particiilar  our  contributions  in  this  paper  are: 

1.  We  continue  the  program  started  in  [3]  by  identifying  a  set  of  conditions  on  the  set  function 
f{S),  for  which  the  parsimonious  property  holds.  In  this  way  we  prove  that  a  large  collection 
of  classical  combinatorial  problems  satisfy  it  including  the  matching  problem,  the  T-joiii 
problem,  a  relaxation  of  the  vehicle  routing  problem,  some  disjoint  path  problems,  some  b- 
matching  problems,  etc.  In  particular  all  problems  considered  in  Goemans  and  Williamson 
[4]  satisfy  it.  We  also  find  that  if  the  set  function  f{S)  does  not  satisfy  this  set  of  conditions, 
the  property  does  not  hold.  We  offer  two  proofs  of  the  property:  a  combinatorial  proof  based 
on  splitting  techniques  originated  in  Lovasz  [8]  and  used  in  Goemans  and  Bertsimas  [3]  and 
an  algebraic  proof  based  on  hnear  programming  duality.  Goemans  [6]  has  also  independently 
developed  this  generalization  using  the  techniques  in  [3] .  The  duality  proof  reveals  a  further 
generalization  of  the  parsimonious  property  to  integer  programming  programs  as  well  (the 
dual  parsimonious  property). 

2.  We  use  the  parsimonious  property  to  prove  interesting  monotonicity  properties  of  the  LP 
relaxation  for  problems  in  this  class. 

3.  We  use  the  parsimonious  property  to  give  genuinely  simple  proofs  of  the  integrality  of  some 
polyhedra  Pj{D):  the  T-join  problem,  special  cases  of  the  Steiner  tree  problem  including  the 
shortest  path  problem  and  the  shortest  path  tree  problem. 

4.  We  further  extend  the  parsimonious  property  under  more  general  conditions  and  examine  its 
implications  in  the  disjoint  path  problem.  We  find  that  this  extension  is  the  source  for  several 
results  in  this  area  and  provides  a  unifying  framework  to  understand  these  results. 


5.  Wc  ofTcr  a  new  proof  technique  that  utilizes  the  parsimonious  property  to  find  bounds  on 
the  ratio  ^  Agv  .  Our  proof  technique  leads  to  a  new  approximation  algorithm  for  this  class 
of  problems  that  compared  with  the  algorithm  proposed  by  Goemans  and  Williamson  [4]  is 
simpler  to  implement  as  it  does  not  use  reverse  deletions,  but  only  shortest  path  computations. 

6.  We  use  the  parsimonious  property  to  prove  new  approximation  bounds  for  Problem  IPf(D) 
with  D  ^H). 

The  paper  is  structured  as  follows.  In  Section  2  we  introduce  the  properties  of  the  set  function 
/(S)  that  imply  the  parsimonious  property  and  examine  classical  combinatorial  problems  that  can 
be  modelled  in  this  way.  In  Section  3  we  prove  the  parsimonious  property  as  well  as  the  dual 
integral  parsimonious  property.  In  Section  4  we  derive  interesting  monotonicity  properties  of  this 
class  of  problems  as  consequences  of  the  parsimonious  property.  In  Section  5  we  further  extend  the 
parsimonious  property  and  apply  it  to  the  analysis  of  the  disjoint  path  problem.  In  Section  G  we 
examine  applications  of  the  peirsimonious  property  to  the  integrality  of  certain  polyhedra  Pf{D).  In 
Section  7  we  introduce  a  new  proof  technique  to  bound  the  ratio  "z^t^-  This  proof  technique  gives 
rise  to  a  new  approximation  algorithm  for  the  problems  considered  in  Goemans  and  Williamson  [4]. 
We  further  examine  applications  of  the  pzirsimonious  property  to  the  approximability  of  Problem 
IPj{D). 

2     Parsimonious  set  functions 

In  their  study  of  the  approximability  of  problems  in  the  class  I Pf  {</}),  Goemans  and  Williamson  [4] 
(for  the  case  that  f{S)  takes  values  in  {0,  1})  and  Williamson  et.  al.  in  [16]  (for  the  case  thai  f{S) 
takes  values  in  Z+)  introduce  the  following  set  of  conditions  for  the  set  function  f{S). 
Conditions  A  (proper  set  functions): 

1.  /(0)  =  O. 

2.  Symmetry:  f{S)  =  f{V  \  S)  for  all  SCV. 

3.  Propereness:  If  /I  n  /?  =  0,  then  f{A  U  B)  <  max{/(/l),  /(/?)). 


We  next  introduce  the  following  set  of  conditions: 
Conditions  B  (parsimonious  set  functions): 

1.  /(0)  =  0. 

2.  Symmetry:  f{S)  =  f{V\S)  for  all  S  QV. 

3.  Node  Subadditivity  (NS):  If  A  n  {x}  =  0,  then  f{A  U  {x})  <  f{A)  +  fi{x}). 

4.  Quasi-supermodularity  (QS):  For  aU  5,T  C  K,  STl  T  ^^  0 
Either 

/(5)  +  /(T)</(5UT)  +  /(5nT) 

or 

f{S)+f{T)<f{S\T)  +  f{S\T). 

We  also  introduce  the  general  subadditivity  condition: 
Subadditivity:  If  /I  n  B  =  0,  then  f{A  U  B)  <  f{A)  +  f[B). 

The  QS  property  was  also  introduced  in  the  recent  paper  of  Goemans  et.  al.  [5].  who  used  llic 
term  weakly  supermodular. 

Finally  we  introduce  a  third  set  of  conditions  on  the  set  function  /: 
Conditions  C  (weakly  parsimonious  set  functions): 

1.  /(0)  =  O. 

2.  Symmetry:  /(5)  =  f{V  \  S)  for  all  S  QV. 

3.  Weak  Subadditivity  (WS):  If  ^  n  {x]  ^  0,  then  f{A  U  {x})  <  f{A).  We  will  then  say 
that  a;  is  a  weakly  Steiner  vertex. 

4.  2-Quasi-supermodularity  (2-QS):  For  every  three  mutually  crossing  sets  (two  sets  A,  B 
are  crossing  li  A\B,  B\A,  Ar\B  axe  nonempty)  at  least  two  of  them  satisfy  the  QS  property. 

Compared  with  Conditions  B  (parsimonious  set  functions),  weak  subadditivity  is  stronger  than 
node  subadditivity,  while  the  2-QS  property  is  a  relaxation  of  the  QS  property.  In  other  words, 
there  are  set  functions  /  satisfying  one  of  conditions  B  or  C  but  not  the  other. 


We  next  show  that  Conditions  B  are  more  general  than  Conditions  A. 

Proposition  1  [Goemans  et.  al.  [5]]  lyct  f  be  a  symmetric,  set  function  with  /(0)  =  0.  Tlien, 
if  f  is  proper,  it  is  quasisupermodular  and  node  subadditive. 

Proof  :  If  /  is  a  proper  function,  then  clearly  /  is  node  subadditive.  Among  the  terms  f{S  Pi  T), 
f{S  U  T),  f{S\T),  f{T\S),  say  f{S  n  T)  attains  the  minimum.  By  properness,  f{S)  < 
max(/(5  n  T),  f{S  \  T))  =  f{S  \  T),  and  f{T)  <  max(/(5  n  T),  f{T  \  S))  =  f{T  \  S).,  and  so 
/(S)  +  fiT)  <  f{S  \T)  +  f{T  \  S).  The  other  cases  follow  similarly  from  symmetry  of  /.  D 

In  Figure  1  (we  consider  symmetric  set  functions)  we  draw  the  relations  of  the  various  conditions 
we  considered.  As  we  show  in  the  next  sections,  the  parsimonious  property  holds  for  set  functions 
/  satisfying  either  conditions  A,  or  B,  or  C. 

Proper 


weakly 
parsimonious 
(WS+2-QS) 


parsimonious 
(NS+QS) 


Figure  1:  Relations  of  parsimonious,  weakly  parsimonious  aind  proper  set  functions. 


Remarks: 

1.  The  following  simple  observation  is  usually  useful  in  checking  whether  a  given  function  has 
the  QS  property.  We  select  any  node  v  and  check  the  QS  properly  only  for  those  sets  S  and 
T  containing  v.  By  symmetry  of  /,  we  can  extend  the  QS  property  to  all  5  and  T. 


2.  Conditions  B  are  strictly  more  general,  i.e.,  there  are  set  functions  which  are  parsimonious 
but  not  proper. 

3.  The  symmetry  conditions  are  without  loss  of  generality.  If  the  function  /  is  not  symmetric, 
we  can  redefine  the  following  symmetric  set  function:  f{S)  =  f{V\S)  =  max[/(5),  /(V''\5)]. 
Notice  that  /Z,-(0)  =  /Z/(0)  and  Z,-(0)  =  Z;(0),  since  the  optimal  solution  of  /P/(0)  is 
feasible  in  IPA^)  and  vice  versa. 

2.1      Examples  of  problems 

In  Table  1  below  we  review  several  classical  combinatorial  problems  formulated  using  the  cutset 
formulation  /P/(0)  for  /  satisfying  both  Conditions  A  and  B.  In  Section  5  we  show  how  Conditions 
C  naturally  arise  in  the  study  of  the  disjoint  path  problem. 


Problem 

ns) 

Conditions 

Spanning  tree 

1  for  all  5  /  0,  K 

A,B 

Steiner  tree 

1  if  5  n  T  7^  0,  T 

A,B 

Shortest  path 

1,  |5n{s,(}|-l 

A,B 

Generalized  Steiner  tree 

1  if  5nTi  7^0,ri,  i  =  l,...,k 

A,B 

Nonbipartite  matching 

1  if  |5|  =2k+l 

A,B 

T-join 

1  if  |5nT|  =2it  +  l 

A,B 

Network  survivability 

max(i,j)e*(S)  ^ij ,  re>0 

A,B 

fc-connected  graph 

k  for  all  5  /  0,  V^ 

A,B 

Tree  partitioning 

1,  if  |5|  <A;,  \S\>n-k 

A,B 

Point-to-point  connection 

1,  if  i5nC|7.^|5nD| 

A,B 

Table  1:  Problems  formulated  as  /P/(0)  satisfying  Conditions  A  (and  therefore  B). 

We  next  describe  problems  that  are  parsimonious  but  not  proper. 

The  b-matching  problem 

Given  numbers  b{i)  such  that  J2ievK''')  =  2^,  the  problem  can  be  modeled  in  the  form  of  IPf{D) 


with  D=V  and 

I      1       l-?|>2,  E.€sfc(0=   2k+l  , 

\  b{i)  5={i},  V\{i} 

Notice  that  the  function  /  is  not  proper,  because  for  A,B  disjoint  whose  union  is  K  \  {i}  the 
definition  is  violated.  However,  the  function  f{S)  is  QS.  While  /  is  not  subadditive  for  general  sets 
A,  B,  it  is  node  subadditive  if  b{i)  e  {a,a  +  \}.  While  in  general  the  6-Tnatching  problem  does  not 
satisfy  the  parsimonious  property,  it  does  satisfy  it  if  b{i)  €  {a,a  +  I}. 
The  capacitated  tree  problem 

Given  a  graph  G  =  {V  L)  {0},  E),  demands  d,,  i  €  V,  a  depot  0,  costs  Ce,  e  e  E  and  we  would 
like  to  design  a  tree  of  minimum  cost  such  that  each  subtree  from  the  depot  has  demand  at  most 
Q.  The  capacitated  tree  problem  is  a  popular  relaxation  of  the  vehicle  routing  problem.  A  valid 
cutset  formulation  of  the  capacitated  tree  problem  is  of  the  type  IPj{%)  witii 


^ 0€5 

It  is  obvious  that  J{S)  does  not  satisfy  Conditions  A,  since  it  is  not  proper  but  it  satisfies  Conditions 
B:  It  is  clearly  symmetric  and  subadditive  as  we  show  below:  For  S,T  cV  such  that  S  f\T  =  (b 

j^S  U  T)  =  2^'^^^^^'  =  /(5)  +  J{T), 

/((Sn  {0}) uT)  =  2^'^^\(^^^)  '  <  /(5n  {0})  <  /(5n  {o})  +  f{T). 
It  is  also  QS,  since  for  S,T  cV  containing  0, 

/(5  U  T)  +  f{S  nT)  =  f{S)  +  f{T). 

The  traveling  salesman  and  vehicle  routing  problem 

The  traveling  salesman  problem  can  be  modeled  as  IPj{D)  with  D  =  V  and  f{S)  =  2  for  all 
S  C  V.  Interestingly,  the  vehicle  routing  problem  can  be  modeled  in  our  framework  as  follows. 
Given  a  graph  G  =  {V  U  {0},E),  demands  dt,  i  e  V,  &  depot  0,  costs  Ce,  e  €  E  and  vehicles 
of  capacity  Q  we  want  to  find  tours  of  the  vehicles  from  the  depot  of  minimum  cost,  such  that 
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the  demand  in  each  tour  does  not  exceed  capacity.  Notice  that  the  capacitated  tree  problem  is 
a  relaxation  of  the  vehicle  routing  problem.  We  can  strengthen  the  formulation  if  we  write  for 
example  f{S)  =  2[     '^^  '].  While  this  set  function  is  subadditive,  it  is  not  QS. 

3     The  parsimonious  property 

The  cut-set  formulation  introduced  in  the  previous  section  captures  many  of  the  classical  opti- 
mization problems  studied  in  the  literature.  It  is  thus  interesting  and  indeed  surprising  that  the 
parsimonious  property  holds  for  the  LP  relaxations  of  these  problems.  In  the  remainder  of  this 
section,  we  prove  the  parsimonious  property  in  two  ways: 

1 .  The  primal  proof  is  an  extension  of  edge  spUtting  techniques  introduced  in  Lovasz  [8]  and  used 
in  [3]  to  prove  the  parsimonious  property  for  the  survivable  network  design  problem; 

2.  The  dual  proof  uses  linear  programming  duality  and  extends  an  observation  of  Frank  [2]  for  the 
matching  problem  to  the  general  class  of  problems  P/(0). 

3.1     A  primal  proof  of  the  pairsimonious  property 

Let  re  be  a  feasible  solution  in  P/(0)  with  x{v,u),x{v,w)  >  0,  where  u,v  and  w  are  vertices  in  G. 
We  split  V  at  {u,  w}  by  some  A  >  0  in  the  following  way  : 

x{v,u)  *—  x{v,u)  —  A 

x{v,w)  <— x{v,w)  —  A 

x{u,  w)  <—  x{u,  w)  +  A. 

The  splitting  operation  will  preserve  feasibihty  of  x  unless  there  exists  a  set  5  C  V  such  that  v  e 
S,u,w  ^  S  and  x{S{S))  =  f{S).  We  call  such  a  set  5  a  tight  set.  We  denote  by  S  the  complement 
of  S.  We  also  use  the  notation  x{6{A,  B))  —  ^e={i,j},ieA,j£B^e  and  x{6{S))  =  x(6{S,  S)). 
We  need  a  preliminary  lemma  regarding  properties  of  tight  sets. 

Lemma  2  If  S,T  are  tight  sets,  and  f  is  QS,  then 


either  (ij  S\T,T\S  are  tight,  x  {  6  {S  n  T,S  U  T))  =  0 
or  (ii)  S  n  T,S  LI  T  are  tight,  x{6{S \T,T\S))  =  0. 


Proof  :    Since  /  is  QS,  wc  first  consider  the  case  J{T  \  5)  +  /(5  \T)  >  f{S)  +  f{T).   In  this 
instance, 

f{T\S)  +  f{S\T)     <    x{6{T\S))+xi6{S\T)) 


=    x{6{S))  +x{6{T))  +2x{6{SnT,SuT)) 


=    f{S)  +  f{T)  +  2x{6{SnT,SuT). 


Uencex{6iS\T))  =  f{S\T),x{S{T\S))  =  f{T\S)  and  x{6{SnT,SuT))  =  0,  i.e.,  condilion  (i) 
holds.  On  the  other  hand,  the  case  f{S  U  T)  +  /(5  n  T)  >  /(5)  +  /(T)  gives  rise  to  condition  (ii), 
using  an  identical  argument.  □ 

We  can  now  prove  the  central  result  of  this  section. 

Theorem  3  (Parsimonious  Property)  //  the  cost  function  c  satisfies  the  triangle  inequality, 
and  f  is  a  parsimonious  set  function,  then 

Zf{D)  =  Z;(0),    for  all  D. 

Primal  proof:  Let  x  be  an  optimal  solution  in  P/(0),  with  5Zee£2;e  minimal.  Suppose  there  is 
a  u  in  D  that  has  x{6{v))  >  f{v);  Let  u  be  such  that  x{v,u)  >  0.  Let  5  be  a  minimal  tight  sot 
containing  v  but  not  u.  (If  such  set  does  not  exists,  then  we  can  decrease  x{v,  u)  by  some  positive 
A,  while  maintaining  feasibility.  This  contradicts  the  minimallity  of  x.)  If  T  is  another  tight  set 


containing  v  but  not  u,  then  x{6{S  nT,S  UT))  >  x{v,u)  >  0.  Since  /  is  QS,  condition  (ii)  of 
Lemma  2  holds  and  so  5  n  T  is  a  tight  set.  By  the  minimallity  of  5,  S  is  contained  in  T.  So  there 
is  a  vmique  minimal  tight  set  containing  v  but  not  u. 

In  addition,  there  exists  a  u;  in  5  with  x{v,w)  >  0,  else  f{S)  =  x{6{S))  =  x{6{v))  +x{6{S\ 
i^}))  >  f{v)  +  f{S  \  {v)),  violating  the  node  subadditivity  of  /.  We  can  then  split  v  at  {u.  w}  by 
some  positive  A,  where 

A<  ^mm{x{6{S))  -  fiS)  :  V  €  S,u,w  ^  S,xi6{S))  -  f{S)  >  0}. 

Because  of  the  triangle  inequality,  this  operation  yields  another  feasible  optimal  solution  x'  witli 
x'{6{v))  <  x{6{v))  and  x'{6{t))  -  x{6{t))  \U  ^v,t  eV.  This  contradicts  the  minimallity  of  x.  □ 
Remarks: 
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1.  In  the  splitting  process,  if  x,f  are  both  even  and  integral,  we  can  choose  A  to  be  1.  This 
corresponds  to  the  classical  edge-sphtting  notion  and  has  played  an  important  role  in  many 
cormectivity  problems.  We  summarize  this  discussion  in  the  following  corollary. 

Corollary  1  Let  G  he  an  Eulerian  multigraph  and  xg  be  the  incidence  vector  of  G.  Lei  f 
be  an  even,  parsimonious  set  function.  If  xg  is  a  feasible  integral  solution  to  IPf{%)  and 
x{d{v))  >  f{{v})  for  some  v  e  V,  then  there  exists  {u,w}  and  an  edge  splitting  operation 
of  V  at  {u,w},  yielding  a  new  graph  G'  and  a  corresponding  incidence  vector  xg'  that  is  a 
feasible  integral  solution  to  IPj{%). 

This  corollary  generahzes  the  following  result  of  Lovasz  [8]  and  a  refinement  due  to  Goemans 
and  Bertsimas  [3]. 

Corollary  2  Let  G  be  an  Eulerian  multigraph,  rG{i,j)  be  the  maximum  number  of  edge 
disjoint  paths  from  i  to  j  and  v  be  a  vertex  of  G.  Then  there  exists  {u,w}  and  an  edge 
splitting  operation  of  v  at  {u,u}}  (obtaining  a  new  graph  G' )  such  that 

•  rG'{i,j)  =  rciij)  ifij  i-  v. 

•  rG'{vJ)  =  mm{rG{i,j),degG{v)  -  2)  for  j  ^^  v. 

Proof  :  Set  f{S)  =  max{rG{i,j)  :  e  =  {i,j)  G  S{S)}.  Note  that  /  is  a  parsimonious 
function.  The  incidence  vector  xg  of  graph  G  is  a  feasible  integral  solution  in  IPf{9),  with 
x{6{v))  =  degciv)  >  f{v).  Since  G  is  Eulerian,  f{S)  and  x{6{S))  are  even  for  each  5  C  V'. 
Thus  there  exists  an  edge  splitting  operation  of  v.  The  graph  G'  obtained  in  this  way  is  a 
feasible  solution  with  xg'{S{S))  >  f{S).  By  the  max-flow-min-cut  theorem,  the  graph  G'  has 
the  required  cormectivity.  □ 

2.  Notice  that  for  the  parsimonious  property  to  hold  /  needs  to  be  QS  and  node  subadditi\'e. 
The  full  subadditivity  is  not  needed.  The  5-matching  problem,  for  example  with  b{v)  =  k 
or  fc  +  1  is  both  QS  and  node  subadditive,  while  if  b{v)  e  {k,k  +  l,k  +  2}  it  is  not  node 
subadditive. 
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Although  the  condition  that  the  cost  function  c  satisfies  the  triangle  inequality  seems  restrictive, 
we  next  show  that  for  problems  of  the  form  /P/(0),  P/(0),  i.e.,  with  no  degree  constraints,  we  can 
ensure  that  this  condition  is  met  by  the  following  transformation.  Let  (/(u,  f)  denote  the  shortest 
path  between  u  and  v  with  c  as  the  length  function.  Clearly  d  satisfies  the  triangle  inequality.  Let 
/Z^(0)  and  Zy-(0)  denote  the  respective  solution  value  with  (/  as  the  objective  function. 

Theorem  4 

Z}(0)  =  Z;(0);    /Z}(0)  =  IZjifD). 

Proof  :  Let  x  be  an  optimal  solution  to  Pf(0)  (or  IP'A^)).  Consider  an  edge  e  =  {u,v)  such  that 
(/(e)  <  c(e).  Let  P  be  a  shortest  path  (with  respect  to  c)  linking  u  and  v.  Then  P  ^  [e],  and 
c!{g)  —  c{g)  for  each  edge  g  on  P.  If  x{e)  >  0,  we  can  reroute  the  flow  through  the  path  P,  resulting 
in  another  optimal  solution  (since  c'{c)  =  d{P)).  Repeating  this  procedure,  we  have  x{c)  >  0  only 
when  c'(e)  =  c(e),  thus  proving  the  theorem.  D 

3.2     A  dual  proof  of  the  parsimonious  property 

The  dual  of  P/(D)  is  as  follows  : 

DZ/(D)  =  maximize  T.scvyiS)fiS) 

subject  to    T,S:eeS(S)  y{S)  <  c(e),  e  e  E 
y{S)>0,ScV,S  ^D. 

Let  DPf{D)  be  the  dual  polyhedron  and  DZf{D)  denote  the  optimal  objective  value.  To  prove  the 
parsimonious  property  using  a  dual  argument,  we  only  need  to  show  that  among  all  dual  optimal 
solutions  to  DPf{D),  we  can  always  choose  one  with  y{v)  >  0  for  all  v  €  D.  This  solution  is  then 
feasible  to  DP}{^).  Let  T  be  a  collection  of  sets  (subsets  of  V).  We  call  this  family  of  sets  laminar 
if  for  all  /t,  /i  €  T  either  Af^B  =  ^,orAcB,OTB(zA. 

Theorem  5  //  the  cost  function  c  satisfies  the  triangle  inequality,  and  f  is  a  parsimonious  set 
function,  then 

Zj{D)  =  Z;(0),    for  all  D. 
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Dual  proof:  Let  y  be  a  dual  optimal  solution  in  DPf{D).  By  the  QS  property,  we  may  assume 
that  the  set  T  :=  {S  :  y{S)  >  0}  is  laminar,  since  we  can  always  replace  two  intersecting  sets  5 
and  Thy  S\T,T\S  or  S  n  T,S  U  T.  Suppose  there  exists  &v  e  D  such  that  y{v)  <  0.  For  all 
A  e  T  containing  v,  we  replace  A  by  V  \A,  i.e.,  we  set 

y{A)  ^  0,  y{V  \  A)  -  y{V  \  A)  +  y{A). 

In  this  way  we  obtain  another  dual  optimal  solution  with  no  member  of  J^  containing  v.  Note  that 
J^  is  still  laminar. 

Let  p(e)  =  Z!s:ee<5(s)  y('S')-  By  dual  feasibility,  p(e)  <  c(e).  We  may  assume  that  there  is  a 
u  £  V  such  that  p{u,v)  =  c{u,v),  since  we  can  increase  y{v)  otherwise.  Let  A  he  a  maximal 
member  of  J^  containing  u.  Let  A  =  imn{—y{v),y{A)).  We  modify  the  dual  solution  as  follows  : 

y{v)    <-    y{v)  +  A 


y{Au{v})     *-    y{Au{v})  +  A 
y{A)    ^    yiA)-A. 

To  check  for  feasibihty  of  this  modified  solution,  we  only  need  to  consider  edges  of  the  form  {v,  w) 
where  w  is  not  in  A.  Note  that  by  the  construction  of  ^,  p{u,w)  —  p{v,w)  +p{u,v)  -2y{v).  Hence 
c{v,  w)  >  c(u,  w)  —  c{u,  v)  =  c{u,  w)  —  p{u,  v)  >  p{u,  w)  -  p{u,  v)  —  p{v,  w)  -  2y{v)  >  p{v.  w)  +  2A. 
Thus  the  modified  solution  is  dual  feasible.  By  repeating  this  procedure,  we  can  construct  a  dual 
optimal  solution  with  y{v)  >  0  for  all  v  in  D.  □ 

Notice  that  if  y  in  the  above  proof  takes  only  integral  values,  then  A  can  be  chosen  to  be  integral. 
This  yields  an  integral  analogue  of  the  parsimonious  property  in  a  dual  sense. 
Let  DIZf{D)  denote  the  optimal  objective  value  over  DPf{D)  with  integrality  constraints  on  y{S). 

Theorem  6  (Dual  Integral  Parsimonious  Property)  ///  is  parsimonious,  andc  satisfies  the 
triangle  inequality,  then  DIZf{il>)  —  DIZf{D). 

3.3     On  the  minimallity  of  conditions  for  the  parsimonious  property 

We  remark  in  this  subsection  that  the  parsimonious  property  does  not  hold  if  we  relax  either  the 
QS  or  the  node  subadditivity  property. 
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Consider  the  set  function  /  on  3  nodes  as  follows:  f{vi)  =  f{{v2,V3})  =  1,  /(5)  =  0  otherwise. 
Then  clearly  /  is  QS,  but  it  is  not  node  subadditive.  In  this  case  the  parsimonious  property  does 
not  hold,  as  the  polyhedron  Pf{V)  is  empty. 

On  the  other  hand,  subadditivity  alone  does  not  guarantee  the  parsimonious  property.  Define  / 
on  4  nodes  as  follows:  f{S)  =  1  if  |5|  =  1  or  3,  f{S)  =  2  otherwise.  Then  /  is  clearly  subadditive 
and  symmetric.  In  this  instance,  Pf{V)  is  again  empty,  since  ii  x{Vi,Vj)  >  0  and  x{vi)  =  x{vj)  =  L 
then  x{6{{vi,Vj}))  <  2. 

4      Monotonicity  properties 

Let  /  be  a  parsimonious  set  function  defined  on  V,  and  let  W  be  a  set  disjoint  from  I'.  Lot 
fw{S)  =  /(5n  V)  for  S  C  V  UW.  We  call  fw  an  extension  of  /  to  V  u  W .  It  is  easy  to  check 
that  fw  is  again  a  parsimonious  set  function.  Note  that  fw{S)  =  0  if  5  C  W .  For  this  reason,  we 
call  \V  the  set  of  Sleincr  vertices. 

For  instance,  when  f{S)  =  1  for  all  S  C  V,  the  extension  fw  corresponds  to  tlie  Steincr  tree 
problem,  with  W  the  set  of  Steiner  vertices.  When  f{S)  =  1  for  all  odd  5  in  V ,  fw  is  the  I '-join 
function.  Goemans  and  Bertsimas  [3]  for  the  survivable  network  design  problem  and  independently 
Shmoys  and  Williamson  [15]  for  the  Held  and  Karp  bound  proved  the  following  monotonicity  result. 
We  extend  this  result  to  the  class  of  parsimonious  functions. 

Theorem  7  Lei  f,g  be  parsimonious  functions  defined  on  V  and  VuW  respectively  and  fw  defined 
as  above.  Suppose  fw{S)  <  g{S)  for  all  S  C  V  L)  W .  If  the  cost  function  c  (defined  on  V  U  \VJ 
satisfies  the  triangle  inequality,  then 

Zf{V)<Zg{VUW). 


Proof 


Zf{V)    =  Zf^iVuW) 

=  Zj^{^)     (parsimonious  property) 

<  Zg(fh)     {g  dominates  fw) 

=  Zg{V  U  W)    (parsimonious  property). 
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n 
Notice  that  the  monotonicity  result  does  not  hold  for  IZf{D)  in  general.  This  is  due  in  part  to 
the  fact  that  the  parsimonious  property  does  not  hold  for  integral  solutions.  On  the  other  hand. 
since  the  parsimonious  property  holds  for  dual  integral  solutions,  we  show  next  that  the  following 
monotonicity  result. 

Theorem  8  Let  f,fw,9  be  defined  as  above  and  fw  <  9-  Jf  c  (defined  on  V  U  W)  satisfies  the 
triangle  inequality,  then 

DlZji^)  =  DlZf^iUi)  <  DIZg{</)). 

Proof  :  Clearly  DIZf^{9)  <  DIZg{ID),  since  fw  <  9-  We  show  next  that  DIZ}{%)  =  DIZj,„{^). 
Let  {y{S)  :  5  C  F}  be  an  optimal  solution  to  DIP}{^).  Let  M  be  a  large  positive  number.  \\'g 
construct  a  feasible  solution  to  DIPj^{W)  as  follows: 

y'{S)  -  2/(5)  if  5  CK, 
y'{v)  =  -M  \{veW, 
y'{S)    =    0  otherwise  . 

Since  fw{v)  =  0  iorv  in  W,  Escvuwy'{S)fw{S)  =  Escv  y(5)/(5).  Therefore,  DIZj^{W)  > 
DIZf{(D).  From  the  dual  parsimonious  property,  DIZf^{<J})  =  DIZf^{W)  >  DIZf{(/)). 

On  the  other  hand,  if  {y{S)  :  S  C  V  U  W}  is  optimal  for  DIZf^{(l>),  then  by  constructing 
y'{S)  ■=  J2T:Tnv=s  2/(^))  we  obtain  a  feasible  solution  to  DIPj{^),  with  the  same  objective  solution. 
Hence  £»/Z/(0)  >  DIZj^{<Ji).  We  conclude  that  DIZ0)  =  DlZf^id).  □ 

5     Weakly  parsimonious  functions  and  the  disjoint  path  problem 

A  natural  question  is  whether  node  subadditivity  and  quasi supermodularity  (QS)  are  the  most 
general  conditions  on  the  set  function  /  for  the  parsimonious  property  to  hold.  In  this  section  we 
show  that  the  parsimonious  property  still  holds  for  weakly  parsimonious  set  functions  (Conditions  C 
introduced  in  Section  2)  and  observe  that  these  relaxed  conditions  provide  a  unifying  understanding 
of  several  results  on  the  the  disjoint  path  problem. 

We  next  prove  that  the  parsimonious  property  holds  for  weakly  parsimonious  functions. 
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Theorem  9  Ia:1  D  be  the  set  of  weakly  Sieiner  vertices.  If  c  satisfies  the  triangle  inequality,  and 
f  is  weakly  parsimonious,  then 

zj{D)  =  zjm. 

Proof  :  The  proof  is  similar  to  the  proof  of  Theorem  3.  Let  v  e  D,u  e  V,  and  suppose  x(6{v))  > 
f{v),x{v,u)  >  0,  where  x  is  an  optimal  solution  in  Pf{9).  Consider  the  minimal  tight  sets  that 
contains  v  but  not  u.  By  the  2-QS  property,  there  exist  at  most  2  such  minimal  sets,  say  S\  and 
52-  Then  all  tight  sets  containing  v  but  not  u  must  contain  one  of  these  two  sets. 

We  show  next  that  there  is  a  it;  in  5i  Pi  Sa  with  x{v,w)  >  0.  Assimiing  the  contrary,  then 

/(S.)  =  x{6{S,))  =  x{S{S,  \  {v}))  +  x{6{{v},S:))  -  xi6{{v},  5.))  > 

f{SA{v})+x{6{{v},S;))-x{6{{v},S,)). 

From  weak  subadditive,  f{Si  \  {v})  >  /(S,);  hence, 

x{6{{v},Sl))<x{6{{v},S,)). 
Since  we  have  assumed  that  x{6{{v},  S\  DS2))  =  0,  we  rewrite  the  inequality  for  z  =  1. 2  and  obtain 


x{6{{v},S2  \  5i))  +  x{6{{v},Si  U  52))  <  x(<5({v},  5,  \  52)). 
and 


xi6{{v},  5,  \  52))  +  xi6{{v},  5,  U  52))  <  x{6{{v},S2  \S,)). 


Hence  x{6{{v},S\  U  52))  <  0,  which  is  a  contradiction  since  x{v,  u)  >  0  and  u  €  5i  U  52.  Therefore, 
there  exists  a  u;  in  S\  n52  with  x{v,  w)  >  0.  By  splitting  at  v  using  u,  w,  we  obtain  another  feasible 
optimal  (because  of  the  triangle  inequality)  solution.  By  repeating  this  procedure,  we  obtain  an 
optimal  solution  in  Pf{D),  thus  proving  the  theorem.  □ 

Similar  to  Corollary  1  the  above  proof  actually  yields  the  following: 

Corollary  3  Let  G  he  an  Eulerian  multigraph  and  Xc  be  the  incidence  vector  of  G.  Let  f  be  an 
even,  2-QS  set  function.  If  xg  is  a  feasible  integral  solution  to  IPj{9)  and  x{6{v))  >  /({;})  /"'" 
some  weakly  Steincr  vertex  v  e  V,  then  there  exists  {u,w}  and  an  edge  splitting  operation  of  v  at 
{u,w},  yielding  a  new  Eulerian  graph  G'  and  a  corresponding  incidence  vector  xq'  thai  is  a  feasible 
integral  solution  to  IPf{(/i). 
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As  we  show  next,  this  corollary  provides  a  unifying  way  to  understand  several  seemingly  unrelated 
results  for  the  edge-disjoint-path  (EDP)  problem. 

5.1      2-QS  functions  and  the  disjoint  path  problem 

Given  an  undirected  graph  G  =  {V,  E),  a  collection  of  source-sink  pairs  {si,  ii}, . . . ,  {sk,  tk},  the 
EDP  problem  asks  whether  there  exists  a  collection  of  edge  disjoint  paths  in  G,  each  joining  a  source 
to  its  corresponding  sink.   Let  H  denote  the  demand  graph,  with  edge  set  {{si,t\), ... ,  {sk,  ik)}- 
Let  xc{e)  =  1  if  e  G  G  and  let  x//(e)  =  1  if  e  G  //. 
Clearly,  a  necessary  condition  for  the  existence  of  these  paths  is  the  cut-criterion: 

XG{S{S))>XHi6{S))   for  all  5  CK. 

There  has  been  an  extensive  literature  (see  for  example  Frank  [2]  and  Schrijver  [13])  that  finds 
conditions  on  G  and  H,  so  that  the  cut-criterion  is  both  necessary  and  sufficient  for  the  existence 
of  a  solution  to  the  EDP  problem.  Let  Kn,Cn  denote  respectively  the  complete  graph  and  the 
cycle  on  n  nodes.  We  also  denote  the  disjoint  union  of  m  copies  of  Kn  by  mKn-  The  following 
results  are  known: 

Theorem  10  If  G  +  H  is  Eulerian,  and  H  is  either  a  double  star  or  a  K4  or  a  C5,  then  tlic 
cut-criterion  is  necessary  and  sufficient  for  the  solvability  of  the  EDP  problem. 

The  case  in  which  //  is  a  2K2  was  proved  by  Rothschild  and  Whinston  [12].  The  double  star 
case  follows  easily  from  their  result.  The  K4  case  was  proved  by  Seymour  [14]  and  Lomonoso\-  [9] 
independently.  The  C5  case  is  due  to  Lomonosov  [9].  See  [2]  and  [13]  for  nice  proofs  and  exposition 
of  these  results. 

At  first  sight  these  results  appear  to  be  unrelated  without  a  unifying  characteristic.  We  could 
use  the  theory  developed  in  this  section  to  identify  the  unifying  characteristic  of  all  the  above  results 
contained  in  Theorem  10.  The  central  reason  is  that  the  set  function  xh{S{S))  in  these  cases  has 
the  2-QS  property.  In  particiilar  it  is  easy  to  prove  the  following  proposition: 
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Proposition  2  The  set  function  xh  (6(5))  has  the  2- QS  property  if  and  only  if  II  does  not  contain 
a  3A'2  or  disjoint  copies  of  K3  and  K^-  This  in  turn  holds  if  and  only  if  II  is  a  double  star  or  a 
K4  or  a  Cs- 

In  order  to  sec  how  Proposition  2  can  be  used  to  prove  Theorem  10,  let  us  rewrite  the  cut 
condition  as  follows: 

XG+Hi6iS))>f{S)  =  2xH{6iS)). 

Under  the  assumptions  of  Theorem  10  and  using  Proposition  2,  xg+h  corresponds  to  a  Eulerian 
graph  (by  assumption),  while  /  is  an  even,  2-QS  set  function.  Let  D  =  V\{s\,ti, . . .  ,Sk,tk}  he  the 
nodes  in  G  that  do  not  belong  to  a  source-sink  pair.  In  our  context  D  is  the  set  of  weakly  Steincr 
vertices.  Applying  Corollary  3  we  can  then  perform  edge-splitting  operations  on  the  edges  of  G  to 
obtain  a  new  graph  G'  that  satisfies  the  cut  criterion,  but  with  edges  incident  only  to  the  sources 
or  sinks.  The  rest  of  the  proof  involves  showing  that  G'  has  the  set  of  edge-disjoiiil-paths  joining 
each  source-sink  pair.  This  follows  from  a  tedious  case  by  case  analysis  which  we  omit  here,  as  it 
is  unrelated  to  the  theme  of  the  paper.  By  reversing  the  edge-splitting  operations,  we  obtain  a  set 
of  edge-disjoint-paths  in  G  that  meets  the  cut  criterion,  thus  proving  Theorem  10. 

6     Applications  in  proofs  of  integrality  of  polyhedra 

An  important  direction  of  research  in  integer  programming  is  the  development  of  techniques  to 
show  integrality  of  the  associated  polyhedra  for  integer  programming  problems.  Perhaps  the  most 
common  proof  technique  is  algorithmic.  Researchers  develop  an  optimal  algorithm  for  a  combina- 
torial optimization  problem,  which  at  the  same  time  shows  integrality  of  a  proposed  formulation 
for  the  problem.  In  this  section  we  show  that  the  parsimonious  property  loads  to  non-algorithmic, 
genuinely  simple  proofs  of  integrality  of  some  polyhedra  Pf{D),  yielding  new  simple  proofs  of  some 
classical  results  as  well  as  some  new  results. 

A  milestone  in  combinatorial  optimization  is  the  proof  of  integrality  (Edmonds  [1])  of  the  perfect 
matching  polyhedron.  This  result  follows  directly  from  the  integrality  of  the  T-join  polyhedron, 
as  the  perfect  matching  polyhedron  is  a  face  on  the  T-join  polyhedron.     Surprisingly,  wo  can 
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derive  the  integrality  of  T-join  polyhedron  from  that  of  the  perfect  matching  polyhedron,  using  the 
parsimonious  property. 

Theorem  11  Let  f{S)  =  1  iJ\Sr^T\  is  odd,  0  otherwise.   Then 

Conv{IPj{^))  =  P/(0). 

Proof  :  Let  fr  be  the  restriction  of  /  to  T,  defined  on  5  C  T.  Note  that  IPf^{T)  is  just  the  perfect 
matching  polyhedron  on  T.  We  show  next  that  IZj{^)  =  2/(0)  for  all  integral  cost  functions  c. 
By  Theorem  4,  we  may  assume  c  satisfies  the  triangle  inequality.  The  following  inequalities  are 
immediate: 

/Z;(0)  <  /Z/^(0)  <  /Z/^(r). 

From  the  integrality  of  the  perfect  matching  problem  IZfj.{T)  =  Zfj,{T).  Prom  the  parsimonious 
property  Zf^{T)  =  Zf^{^)  =  Zf  {(!)),  yielding  that 

/Z/(0)  <  Zf{<D). 

The  reverse  inequality  holds  trivially  and  so  IZj{^)  =  Zj{^),  which  shows  integrality  of  the  T-join 
polyhedron.  □ 

In  the  next  section  we  briefly  review  another  (and  in  our  opinion  quite  powerful)  proof  technique 
that  proves  the  integrality  of  the  perfect  matching  polyhedron  directly. 

The  shortest  path  polyhedron  can  be  treated  as  a  Steiner-1-cormectivity  polyhedron  on  two 
terminal  nodes.  Integrality  of  the  polyhedron  also  follows  easily  from  the  parsimonious  property. 
We  generalize  this  result,  using  the  parsimonious  property,  and  show  that  the  cut  set  formulation 
for  the  Steiner-2-Connected  polyhedron  with  at  most  5  terminal  vertices  is  integral. 

Theorem  12  For  the  Steiner-2-Connected  problem  on  at  most  5  terminal  nodes, 

Conv{IP0))  =  PfiH)). 

Proof  :  It  is  well  known  ([11])  that  the  TSP  polyhedron  on  at  most  5  nodes  is  integral.  From  the 
parsimonious  property  the  result  follows  easily.  D 
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We  next  consider  the  multi-commodity  flow  problem  with  a  single  source  s,  multiple  sinks  {t\.t2-  ■  ■  ■ ,  tk} 
and  with  no  capacity  constraints.  Let  D  =  {s,t\,t2,  ■  ■ .  ,tk}-  Algorithmically,  the  problem  reduces 
to  the  computation  of  the  corresponding  shortest  paths  between  the  source  and  the  sinks.  Note 
that  if  D  =  V^  the  problem  is  the  shortest  path  tree  problem.  We  show  that  the  cut-set  formulation 
for  the  problem  with  f{S)  =  Y.u&s  1  if  5  ^  5  and  /(5)  =  f{S)  is  integral.  This  result  also  follows 
from  Johnson  [7]. 

Theorem  13  For  the  uncapaciiated  multi- commodity  flow  problem  with  a  single  source  and  mul- 
tiple sinks, 

C(mv{IPf{<ll))  =  P}{<Ji). 

Proof  :  We  only  need  to  show  that  I  Zj{D)  =  Zf{D)  when  c  satisfies  the  triangle  inequality.  Since 
PfiY)  has  only  a  single  integral  solution  with  x{s,  t,)  =  1  for  each  z  =  1, 2, . . . ,  A:,  the  result  follows 
immediately.  D 

7     Applications  in  worst  case  analysis 

In  recent  years  there  has  been  a  lot  of  interest  in  the  approximability  of  combinatorial  optimization 
problems.  Typically  researchers  propose  a  heuristic  algorithm  for  an  integer  programming  problem 
(a  minimization  problem)  and  compare  the  value  of  the  heuristic  to  the  value  of  the  LP  relaxation 
(or  to  the  value  of  a  dual  feasible  solution  of  the  LP  relgixation).  A  very  nice  and  \ory  general 
example  of  this  approach  is  the  2(1  -  ^)  approximation  algorithm  (T  =  [v  6  V^  :  f{v)  —  1}) 
proposed  in  Goemans  and  Williamson  [4]  for  the  problem  I Pf{9)  with  /  being  proper  (Conditions 
.4)  and  taking  values  in  {0,1}.  A  corollary  of  their  result  is  the  bound  -^^v  ^  2(1  -  -rfi).  A 
distinct  characteristic  of  their  method  is  a  reverse  deletion  step,  in  which  edges  that  were  added 
in  the  solution  are  deleted.  Moreover,  for  the  matching  problem  the  bound  is  exact  (the  matching 
polyhedron  Pf{V)  is  integral). 

In  this  section  we  propose  a  new  proof  method  that  shows  that  ^  ■^L  <  2(1  -  4^)  for  proper 
functions.  The  proof  method  gives  rise  to  a  new  (and  in  our  opinion  more  natural)  algorithiii  thai 
does  not  use  reverse  deletions  and  therefore  it  is  easier  to  implement.  Moreover,  wo  remark  that 
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the  proof  method  is  quite  powerful  as  it  can  prove  integrality  of  the  matching  polyhedron.  It  can 
also  be  used  to  prove  the  integrality  of  the  multicut  formulation  for  the  minimum  spanning  tree 
problem  and  the  branching  polyhedron.  Finally  we  use  the  parsimonious  property  to  bound  ^  Yu/ 
if  c  satisfies  the  triangle  inequality. 

7.1     A  proof  technique  to  bound  the  ratio  -^4^ 

We  consider  problem  /P/(0)  with  /  being  a  0  —  1  proper  function.  Let  T  =  {v  e  V  :  f{v)  =  1}. 
Our  proof  technique  uses  the  crucial  observation  that  a  minimal  solution  to  the  problem  must  be 
a  forest,  and  thus  has  at  most  IT"!  —  1  edges. 

Theorem  14  (Goemans  and  Williamson  [4])   If  f  is  a  0  —  1  proper  function 

iZf{<ii)<2{i--^^)Zjm. 

Proof :  For  the  purpose  of  contradiction  we  assume  the  contrary.  Therefore,  there  exists  a  coimter- 
example  on  the  least  number  of  nodes,  with  /  proper  and  c  integral.  We  may  further  assume  that 

5Ze€£;'^(^)  's  minimal. 

Suppose  there  is  a  u  with  f{v)  =  0.  Let  /'  denote  the  restriction  of  f  on  V  \  {v} .  It  can  easily 
be  checked  that  /'  is  still  proper.  By  the  minimallity  of  the  counter-example, 

Since  the  optimal  solution  in  /P/'(0)  is  also  feasible  in  /P/(0),  /^/(0)  <  IZf  {</}).  From  Theorem 
4,  by  using  the  shortest  path  distances  Z/(0)  =  Z^(0)  and  Z//(0)  =  Z^,(0).  But,  Z^,(0)  =  Z'^{{v}). 
By  the  parsimonious  property  Z^(0)  =  Zj{{v})  =  Z^/(0).  Therefore, 

/Z;(0)<2(l-|ii)Z;(0), 

which  is  a  contradiction.  So  we  may  assume  f{v)  =  1  for  all  v. 

If  there  is  an  edge  e  =  {u,v)  e  E  with  Ce  =  0,  then  by  contracting  this  edge,  and  treating  {u.  v} 
as  a  supernode,  we  restrict  the  problem  to  one  of  strictly  smaller  size.  By  the  minimallity  of  the 
coimter-example,  there  exists  a  solution  that  satisfies  the  theorem.  By  introducing  the  edge  (u.v)., 
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with  no  extra  cost  since  Ce  =  0,  if  necessary,  wc  obtain  a  solution  feasible  to  the  original  problem 
and  the  theorcnn  holds.  Therefore,  we  may  assume  Ce  >  0  for  all  e. 

Now  let  y{v)  =  ^  for  all  v  and  consider  the  cost  function  c'  where  c^  =  Ce  - 1  (Cg  >  0  from  the  pre- 
vious paragraph).  By  the  mininiallity  of  c,  there  exist  x,  y'  such  that  x  £  JPf{9),  5Zs:e6«(S)  v'i'^)  - 
4  and 

J2cf,xe<2{i-^^)  Y,  y'is)f{s). 

Since  /  is  0  -  1,  le  corresponds  to  a  forest  and  therefore, 

X;  Xe  <  iri  -  1  =  2(1  -  i^)  ^  y{v)  =  2(1  -  ^)  X:  f{v)y{v). 

The  last  equality  holds,  since  we  have  shown  that  wc  can  assume  f{v)  =  1. 
Let  y'  =  y'  +  y.  Note  that 

E     y'{S)=     E     y{S)  +  l  +  \<c',  +  \=Ce 

S:e€S{S)  S:ee6{S) 


and  so  y*  is  dual  feasible.  Therefore, 

\T\ 


^   CeXe  =  5:   c'^Xe  +  E   ^^  <  2(1  -  -)-)  J]    y'{S)f{S). 


This  is  again  a  contradiction  and  the  theorem  follows.  □ 

Remarks: 

1.  The  dual  variables  y  constructed  in  the  proof  arc  half-integral.  We  call  a  cut  6(S)  an  /-cut 
if  f{S)  >  0.  We  can  refine  the  previous  theorem  as  follows.  We  have  shown  that  IZf{9) 
is  bounded  above  by  2(1  -  tL)  times  the  maximtmi  half-integral  c-packing  of  /-cuts.  This 
observation  has  an  interesting  implication  for  the  TSP.  It  is  well-known  that,  if  the  cost 
function  c  satisfies  the  triangle  inequality,  the  Christofides  heuristic  constructs  a  solution 
with  objective  value  (denoted  Zc)  not  more  than  3/2  times  of  the  optimum.  This  ro.'^ult 
has  been  strengthen  further  by  Wolsey  [17]  and  Shmoys  and  Williamson  [15]  who  showed 
that  Zc  <  (^  -  ^)Zy(0),  where  /  corresponds  to  the  TSP  function.  We  can  strengthen  the 
inequality  by  replacing  Z/(0)  with  the  vahie  of  the  maximum  half-integral  c-packing  (denoted 
by  DZ;(l/2)).   Note  that  all  cuts  are  /-cuts  in  this  instance  and  2  DZ/(l/2)  <  DZ;(0)  = 
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Zf{0).  From  the  above  discussion,  the  solution  to  the  minimum  spanning  tree  is  bounded 
above  by  2(1  —  ^)DZf{l/2).  A  well  known  result  on  matching  (see  [10])  says  that  the 
minimiun  matching  on  the  set  of  odd  nodes  is  bounded  above  by  DZj{\/2).  Hence 

^c<(3-||^)DZ;(l/2). 

2.  The  proof  only  works  for  proper  and  not  the  more  general  parsimonious  functions,  because 
we  want  the  parsimonious  property  to  hold  even  if  we  create  supemodes.  Therefore,  node 
subadditivity  is  not  sufficient.  Therefore,  we  need  /  to  satisfy  the  full  subadditivity  and  the 
QS  property,  which  for  the  case  of  0  —  1  functions  is  exactly  the  class  of  proper  functions. 

3.  The  above  proof  technique  can  be  used  to  prove  the  integrality  of  the  matching  polyhedron, 
the  multicut  formulation  for  the  minimimi  spanning  tree  problem  and  the  branching  polyhe- 
dron. For  the  matching  polyhedron  the  difficult  step  is  the  case  with  Ce  =  0,  which  can  be 
handled  using  techniques  from  [10].  The  final  step  is  easy,  since  Yle^e  =  V  ~  ^i'  y(^')- 

As  an  example  of  a  different  application  of  the  proof  method  let  us  consider  the  muUicTit 
formiilation  of  the  MST.  Let  n  =  {S\, . . .  5|n|}  be  a  partition  of  V. 

IZmcut  =    minimize  HeeE  CeXe 

subject  to    T.ees(Si,Sj);  i<t<j<|n|a;e  >  |n|  -  1,  V  n  =  {5i, . . .  5|n|} 

Xe  e  {0,1}, 

Let  Zmcut  be  the  LP  relaxation  and  consider  the  dual  problem. 

Zmcut  =    maximize  End^l  -  l)y(n) 

subject  to    En:  eesiSi,Sj)  vi'^)  <  Ce,  ye  £  E 

2/(n)  >  0, 

We  need  to  show  that  IZmcut  =  Zmcut-  We  use  an  identical  proof  method,  assuming  the 
existence  of  a  minimal  counter-example.  The  only  difference  is  that  we  update  the  dual 
variables  as  follows.  Let  11  ==  {{1},  {2}, . . . ,  {n}}  and  y{U)  =  1.  Note  that  J2e^e  =  n  -  \  = 

E(|n|  -  i)2/(n). 
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AllhouKh  the  proof  method  in  Theorem  7.1  is  non-algorithmic,  it  also  leads  to  an  algorithmic 
method  to  construct  an  approximate  solution.  It  differs  from  the  Goemans  and  Williamson's 
algorithm  in  that  it  noods  a  pre-processing  step  to  compute  pairwise  shortest  paths.  With  this  in 
hand,  we  can  discard  all  vertices  with  f{v)  =  0.  We  call  these  vertices  the  Steiner  vertices.  This 
approach  avoids  the  critical  reverse  deletion  step  in  the  Goemans  and  Williamson's  algorithm,  at 
the  expense  of  computing  pairwise  shortest  paths.  Our  algorithm  is  as  follows: 
Approximation  Algorithm  for  0—1  proper  functions 

1.  Compute  the  pairwise  shortest  path  distances  for  all  pairs  of  non-Steiner  nodes. 

2.  Discard  the  set  of  Steiner  nodes.  Select  an  edge  with  the  least  cost.  Merge  the  two  end  nodes 
into  a  supcmode,  delete  all  edges  joining  these  two  nodes  and  update  the  costs  of  joining  the 
supemodes. 

3.  Repeat  step  2  until  two  supemodes  merge  to  form  a  set  S  with  f{S)  =  0.  Let  e'  be  the 
last  edge  selected.  If  there  are  no  more  non-Steiner  nodes,  go  to  step  4.  Else  for  all  edges 
remaining,  reduce  the  cost  to  c(e')  and  return  to  step  1. 

4.  Replace  the  edges  selected  in  Step  2  and  3  by  its  corresponding  shortest  paili  in  the  original 
graph  G. 

For  the  Steiner  tree  problem,  our  algorithm  emulates  the  MST  heuristic  on  non-Steiner  nodes 
with  the  pairwise  shortest  distaince  metric.  In  this  respect  Theorem  7.1  generalizes  the  well  known 
fact  that  the  MST  heuristic  gives  a  2(1  -  tjtt)  approximate  solution  to  the  minimum  Steiner  tree 
problem. 

For  arbitrary  proper  functions  /,  as  Goemans  and  Williamson  [4]  observe,  we  can  construct 
a  feasible  solution  by  utilizing  Theorem  7.1.  Let  p\  <  p2  <  . . .  <  Pn  be  the  distinct  values  of 
/,  and  for  each  i,  fp,{S)  =  1  if  f{S)  >  p,  and  0  otherwise.  Note  that  /p,  is  proper  0  -  1.  By 
appending  p,  -  p,_i  {po  =  0)  copies  of  the  approximate  solution  to  /p,  for  each  i  =  1, 2, ....  n.  we 
obtain  a  feasible  solution  which  is  within  2H(pi,p2, . . .  ,Pn)  times  of  the  optimal  solution,  whore 
'W(Pi.P2)  ■■■  ,Pn)  =  J17=\  '''"J''"'  ■  Similarly  for  arbitrary  QS  functions  we  can  use  llic  results  of  [5] 
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to  find 

/Z;(0)<2W(p,,P2,...,Pn)Z/(0). 

7.2     On  the  approximability  of  IPf{V) 

We  next  study  the  approximability  of  IPf{V).  The  recent  research  activity  on  approximation 
algorithms  has  so  far  concentrated  on  problem  /P/(0),  partly  because  it  is  difficult  to  construct  a 
feasible  integer  solution  to  IP/ {V).  In  fact,  checking  feasibility  is  usually  NP-hard,  as  indicated  for 
the  case  of  the  Hamiltonian-Cycle  problem.  Using  our  understanding  of  edge-splitting  techniques 
and  the  parsimonious  property,  we  can  extend  many  of  the  approximation  results  to  IZf{\').  when 
/  is  an  even  parsimonious  function,  and  c  satisfies  the  triangle  inequality. 

Theorem  15  ///  is  an  even  parsimonious  function,  and  c  satisfies  triangle  inequality,  then 

IZf{V)<2n{puP2,...,Pn)Zf{V). 

Proof  :  Let  /'  be  f/2.  Then  /'  is  again  a  parsimonious  function.  From  [5]  We  first  construct  an 
approximate  solution  to  JPf'{9),  with  x\y  denoting  the  primal  and  dual  solution  respectivel}'.  Let 
X  =  2x'.  Then  the  graph  corresponding  to  x  is  Eulerian,  since  each  vertex  has  even  degree.  Note 
that 

n{^,...,^)=n{p,,...,pn),  and 
^  CeXe  =  2Y^  c,x',<An{^,...,^-)Y.  t/(5)/'(S)  -  27f  (p, , . . . ,  p„)  ^  y{S)f{S). 

e  e  ^  ^        S  S 

Since  /  is  even  and  the  graph  corresponding  to  x  is  Eulerian,  applying  Corollary  1,  we  can  use 
edge-splitting  operations  to  construct  a  feasible  solution  to  IZj{V).  By  the  triangle  inequality,  the 
cost  of  the  constructed  solution  has  not  increased  and  therefore, 

/Z/(K)<27i(pi,P2,...,p„)Z;(0). 

From  the  parsimonious  property  Z^(0)  =  Zf{V)  and  the  theorem  follows.  □ 

Remark:  For  the  case  of  the  TSP,  the  previous  theorem  corresponds  to  the  well-known  fact  that 

doubling  the  edges  of  the  MST  solution  yields  a  2-approximate  solution  to  the  TSP. 
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